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ABSTRACT 

The equation of motion for an elastically mounted, linearly damped 
circular cylinder subjected to a time-dependent transverse fluid force 
has been formulated. The steady-state and transient solutions have 
been obtained in terms of the in-phase and out-of-phase force coeffi- 
cients, amplitude ratio, and the internal damping of the cylinder. The 
results have shown that a cylinder may be self-excited in the range of 
Strouha I frequencies close to the natural frequency of the cylinder. 

It has further been shown that the heuristic models previously used are 
inadequate to correctly predict the in-phase or the inertia component 
of the exciting force. Additional experiments and analysis could en- 
able one to predict the behavior of cables subjected to ocean currents. 
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NOMENCLATURE 



A 

C 



m I 
D 

F 

T 

t 

U 

U 

m 

V 

A/D 

D/VT 

VT/D 

Re 

S 

f 

s 

f 

c 

V 

Pf 

P. 



amp I i tude of cy I i nder osci I I at ion 
drag coefficient as defined by equation (29) 
inertia coefficient as defined by equation (30) 
Fourier-averaged drag coefficient 
Fourier-averaged inertia coefficient 
diameter of the cylinder 

instantaneous total force acting on the cylinder 

period of oscillation 

time 

instantaneous velocity 

maximum velocity in a cycle (U^ = 2 tt A/T) 

mean velocity of the uniform flow 

dimensionless amplitude ratio 

dimensionless frequency parameter 

dimensionless period parameter 

Reynolds number (Re = VD/V) 

Strouhal number (S = f D/V) 

s 

vortex-shedding frequency 
cy I i nder osci I lation frequency 
phase ang I e 

fluid kinematic viscosity 
fluid density 
cy I i nder dens i ty 
density ratio, P^/Pg 



7 



ACKNOWLEDGEMENT 



I would like to express my great appreciation to Professor T. 
Sarpkaya for his guidance throughout the period of research. 

It is also a pleasure to acknowledge the help given by the per- 
sonnel of the machine shop of the Department of Mechanical Engineer- 
ing in constructing the testing apparatus. In particular, I wish 
to acknowledge the efforts of Mr. J. McKay in the construction of 
the equipment. 



8 



I. INTRODUCTION 



Experience has shown that elastic structures of one or more degrees 
of freedom, such as tall slender structures, cables, etc., can extract 
energy from the flow about them and can develop catastrophic flow- induced 
oscillations. The understanding of this energy-extraction process is of 
paramount importance if one is either to eliminate or minimize it or to 
design the elastic structure such that it can withstand the oscillations 
under the contemplated environmental conditions. 

Evidently, the real origin of the self-excited oscillations lies in 
the basic mechanism of the vortex-shedding process. Consequently, the 
prediction of such oscillations requ i res a tractable theoret i ca I -formu I a- 
tion of the vortex-shedding phenomenon. The fact that such a formulation 
does not exist even for separated flow about a stationary cylinder shows 
the complexity of the treatment of flow about oscillating bodies. The 
problem would have been relatively simpler had the vortex shedding and 
the forces arising from it been independent of the motion of the body. 
Experiments show that vortex generated forces give rise to oscillations, 
and the oscillations modify the forces, and so on, finally resulting in 
a steady-state oscillation such that the Strouha I frequency of the vortex 
shedding nearly coincides with the natural frequency of oscillation of 
the body. 

Mathematical models of flow-induced vibrations of bluff bodies and 
the response of circular cylinders to vortex shedding have been aptly 
described by Parkinson Clj, and Currie, et a I . C2] and will not be 
repeated here. 
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It appears that among the various models considered so far, the 
"wake-oscillator" model of Hartlen and Currie C3H attracted considerable 
attention among the students of vibration analysis. This model predicts 
some of the observed features of the self-excited oscillations even 
though the reasons as to why an osci I lator of the van der Pol type 
should at all describe the self-excited oscillations of an elastic body 
are not yet clear. There are several coefficients in the model which 
must be determined experimentally. There are, at present, wide varia- 
tions in the methods of measurement of some of the coefficients. For ex- 
ample, the damping of the vibrations of a given beam or cylinder is 
determined in still air or water depending on whether the self-excited 
oscillations take place in air or water. The added mass is assumed to 
be equal to the displaced mass, of the cylinder. There is no distinction 
between the material damping as would be determined in vacuum, and the 
damping caused by the fluid. 

It became clear, during the course of the present investigation, 
that the forces acting on a cylinder oscillating harmonically transverse 
to a uniform flow should be measured and the in-phase and out-of-phase 
components of the force should be determined. This phase of the inves- 
tigation has been carried by Meyers C4H and Fortik CSH under the 
direction of Professor T. Sarpkaya, and the drag and inertia coeffi- 
cients for the lift force have been determined through the use of 
Mori son’s C63 equation. 

The purpose of the work reported herein is to make use of the 
measured force by combining it with the equation of motion of an 
elastically mounted, linearly damped, circular cylinder and to attempt 
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to predict the self-excited oscillations of the cylinder for given and 
directly measurable physical parameters such as the mass, length, 
diameter, material damping, and the natural frequency of the body, and 
the density and velocity of the fluid. 

It will be assumed that the measurements and the predictions are 
in the same Reynolds number range and that the force coefficients do 
not measurably depend on the Reynolds number in the range under considera- 
tion. 
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II. ANALYSIS 



A. FORCING FUNCTION 

The transverse force acting on a circular cylinder undergoing har- 
monic oscillations in the transverse direction to a uniform flow is 
expressed in terms of the Morison’s equation as C4H 



C, = 






^ 2,m.,D.2 . 277, 

C I TT (— tT-) ( — ) sin ^ t 

ml D T 

_ /^m'''.2 . D .2 I 2tt , I 2tt , 

- C ,, (— p— ) ( — ) cos -=- t cos t 

d I U ' I ' I 



( I ) 



where C^| and C^ | are the inertia and drag coefficients, respectively. 
Evidently, the normalized force as we i I as the drag and inertia co- 
efficients depend on the parameters U^T/D = 2ttA/D, and D/VT. The first 
parameter is directly proportional to the relative amplitude, and the 
second parameter may be regarded as the frequency parameter (D/VT) or 
the reduced velocity (VT/D). 

For the purpose under consideration it is more advantageous to ex- 
pand the I cos (jot|cos oot term in series and to retain only the first 
term. This procedure yields, 

^ _ r- 2.^m''',,D.2 . , 8 . ,^m^.2.D,2 ^ 

C, - C . TT (— ( — ) Sin cot - pr— ) ( — ) cos (ot (2) 

L m I D r;-p 3tT d I D r;_ 



Equation (2) could be simplified further by combining it with the equa- 
tion of motion for an elastically mounted cylinder. 
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B. ELASTICALLY MOUNTED CYLINDER 




Fig. I. Elastically mounted, linearly 
damped cylinder. 



The equation of motion for the elastically mounted, linearly damped, 
and periodically forced cylinder may be written as (see Fig. I) 



mx + cx + kx 



F = ^Pf 



2 

DL V C|_ 



( 3 ) 



where m represents the mass of the cylinder, c the linear damping, and 

k the spring constant. The derivatives of x are taken with respect to 

real time t, i.e., — x, etc. 

dt 

Def i n i ng. 



OJ 



= /k/m , 



(1) = 



Q = 



2m ( 1 ) ' o 



f /f , a 
c n ' 



2tt f , t = (0 t , 
n n 



f /f = SV/f D , 
s n n ' 

2 2 
, LD^ f 
pf s 

„ 2 - 22 ' 

2m a; S oj 
n o 



x/D = x , 
r ' 

X = A/D 
r 
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equation ( 3 ) may be reduced to 



? ? OS 7 

X + 2?x + X = ao3 [C , 7T (2 ttX )(^)^ sin Qx 

^ ^ ^ omi roj 



r r 



A c (2ttX (— )^ cos nx] 

3ir d I r (jj 
o 



( 4 ) 



3 2 

Defining noting that = 2ir aS , equation (4) re- 

duces to 



7 I 

X + 2^x + x = pXfi(C.sinnx 7 X C , . cos fix) 

r r r r ml , 2 r dl 

3tt 



(5) 



or by a change of variable from to = x/A, and further simplifica- 
tion one has 



y + 2<;y -I- y = p n^(C . sin fix ^ X C , , cos fix) 

V 'r V ^r ml , 2 rdl 

3tt 



( 6 ) 



Equation ( 6 ) may be written as 



2 

y -t- 2i;y -t- y = Efl sin (flx - a) 

; r p ' p 



(7) 



in which 



/ 






( 8 ) 



and 



I 6 A C 



a = arc tan 



dl 



3tt^ D C 



ml 



(9) 



The coefficient Efl may also be expressed in terms of the maximum 
lift coefficient defined by 



F = i P, V^DL (C, ) sin (fix - a) 
2 f L m 



( 10 ) 
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Inserting equation (10) into equation (3) and simplifying one has 



EQ 



2 



7 

2 

m A 0 ) 



(II) 



or 



E = P. 



VT 2 

2.^ 



(C, ) 

L m 

(A/D) 



( 12 ) 



Combining equations (12) and (8) and simplifying one has 

<=L'm 21 , 

TaTdT- ^ 

It should be noted that if (A/D) j were independent of A/D, then 

it is easy to show that (C, ) /(A/D) would also be independent of A/D 

since the experiments have shown that C^| is independent of A/D. This 

would make the results independent of A/D, a conclusion which is contrary 

to observations. This proves that (A/D) C_,, and (C, ) /(A/D) are indeed 

d I L m 

dependent on A/D in a complex manner. But this dependence is rather 
weak. That is why slight changes in, for instance, can lead to large 
changes in or 

C. SOLUTION OF THE EQUATION OF MOTION 

As shown above, a convenient form of the equation of motion is 

y + 2Cy + y = sin (f^T - a) (14) 

' r ' r ' r 

I . Transient Solution 

The Homogeneous equation correspond i ng to equation (14) is 

+ 2Cy^ + = 0 (15) 



)■ 






^ml 



+ (- 



3tt 



D ^dr 



(13) 
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Its characteristic equation reduces to 



s^ + 2cs + I = 0 (16) 

whose roots are given by 

S| 2 = -^ ± i/H^ (17) 

Considering only the underdamped case, it follows that ? < I. 
Therefore, according to Boyce and Oi Prima C73, the solution to the 
correspond i ng homogeneous equation, that is, the transient solution is 
given by 

(y K = Y. sin (/l - T - B) (18) 

'r h h 

where Y, and 3 wi I I be determined from the initial conditions, 
h 

2, Steady-State Solution 

The particular solution to equation (14) is a steady-state 
oscillation of the same frequency of the excitation. Thus, as indicated 
in Thomson Cs3, the particular solution is assumed to be of the form 

(y ) = Y sin [(^x - a) - 4>] (19) 

f p p 

where (J) is the phase of the displacement with respect to the excitation, 
as shown in Fig. 2, and is given by 




Fig. 2. Vector diagram. 
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and 






= arc tan — ^ 






( 20 ) 



and Y is given by 



Y = (21) 

P ± A\-Q^)^ + (2?n)^ 



Adding the transient and the steady-state solutions, the following com- 
plete solution is obtained 



y = Y^ e~^'^ sin (/l-?^ Ttg) + Y sin 
V h p 



(flT-a-(j)) 



( 22 ) 



Also, its derivative with respect to x is given by 



y = Y {-?e sin ( /l x+B) + e /l-?^ x+B)} 
r h 

+ Y fl cos (fix-a-(j)) (23) 

P 

Applying the first initial condition that x = 0 corresponds to y^ = 0, 
one finds 

Y sin (-a-4)) (24) 

Y = - J — 

h s i n B 

Applying the second initial condition that x = 0 corresponds to 
one obtains 

0 = Yj^(-? sin B + cos B) + Y^fl cos(-a-<p) (25) 

Inserting the expression for Y^ given by equation (24) in equation (25), 
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one has 



B = arc tan 



( 



h-r:,^ 

i;+fl cot(-a-<j)) 



) 



Also, reca 1 1 i ng that as T x -»■ -A si n B, and that 

it is seen from equation (22) that = -1. Introducing 
equation (21) and noting that E _> 0, one obtains 

Substituting the expressions for and Y^ into equation 
tains the following complete solution of the equation of 

sin(-a-<j)) -C,x . , T. TT 

y = : — e ^ sin xtB) -sin (f2T-a-(b) 

, r s I n B ■t. K T 

The results obtained through the use of equation 
presented in connection with the discussion of results. 



(26) 

= x/A -sin 9 , 
this value in 

(27) 

(22) , one ob- 
mot i on 

(28) 

(28) will be 
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III. DISCUSSION OF RESULTS 



The time dependent force acting on the cylinder was expressed, as 
indicated previously, in terms of C , and C,,. Assuming this force to 
be applied to an elastically mounted and linearly damped cylinder, the 
equation of motion for the cylinder was developed and solved in terms 
of the given parameters and experimentally obtained coefficients. Two 
of these coefficients, and inertia and drag coefficients respec- 

tively, were introduced in the expression of the time dependent force 
acting on the cylinder. Using the experimental apparatus and procedure 
described in [l4ll, data were obtained, which together with the results 
from F43, were used to plot C , and C,, as a function of VT/D for 
A/D = 0.25, 0.50, 0.75 and 0.84, as shown in Figs. 3 and 4. 

The results plotted in Fig. 3 show that for the region of interest, 
i.e., for VT/D in the range of 4 to 7, and for the precision of the 
experiments, C^| is independent of A/D. 

On the other hand, j is dependent on both VT/D and A/D as shown 

in Fig. 4. As noted previously, (A/D) C , , and (C, ) /(A/D) are not con- 

d I L m 

stants. Furthermore, it can be shown that (C, ) /(A/D) decreases as 

L m 

A/D increases. Figures 5, 6, 7 and 8 show (C^)^ as a function of VT/D 
for A/D = 0.25, 0.50, 0.75 and 0.84 respectively. 

Therefore, an iterative procedure is needed to determine the complete 
solution of the equation of motion for a specific case, as indicated 
below. 

Firstly, for a given VT/D, C^| is obtained from Fig. 3. Then, for 
a specific relative amplitude A/D and the same VT/D, | is obtained from 



19 



Fig. 4, and the value of is obtained from Figs. 5, 6, 7 or 8, or 

alternatively from equation (13). Then, E is calculated using equa- 
tion (12), since p is known. Substituting this value for E in equa- 

’’ VT 

tion (27), a value for ^ is calculated, since n = (— ^)/(-^). A plot 

of I/E versus fl for specific values of C is shown in Fig. 9. Then, 

the values for a, (j) and S are calculated from equations (9), (20) and 

(26), respectively. Substituting these results into the equation of 

motion and plotting the results, one obtains curves similar to those 

shown in Figs. 10, II and 12. These figures show clearly the effect of 

C on the response of the system. A lowly damped system results in an 

initial overshoot of the amplitude before it reaches a steady-state. 

A highly damped system, on the other hand, reaches its steady-state 

gradually from its initial rest position. 

Evidently, one can perform various parametric studies on the effect 
of damping, density ratio, etc., on the transient and steady-states of 
the oscillation provided that sufficient experimental data are avail- 
able to obtain accurately enough the drag and inertia coefficients. It 
is also evident that the practical problem is the determination of 
whether a cylinder of given natural frequency, density ratio, internal 
damping, etc., will be self excited in a fluid of given velocity. 

The determination of this problem requires an iteration not only 
over VT/D but also over A/D in order to obtain the appropriate values 
of the drag and inertia coefficients. Preliminary attempts to resolve 
this aspect of the problem have resulted in the computer program given 
in Appendix A. This program needs to be improved with additional data 
particularly to the lower and higher values of A/D where the force 
transfer coefficients vary most rapidly with A/D. 
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Then, it will be possible to improve the program to predict with 
greater precision the range of self-excited oscillations as well as 
the i r amp I i tude. 

A computer program for the complete solution given by equation (28) 
and for the plotting of the results is given in Appendix B. 

It is evident from the foregoing discussion that the prediction of 
the characteristics of the self-excited oscillations is based on the 
measured time-dependent transverse force rather than on an assumed 
variation of the lift coefficient through the use of a heuristic model 
C3H. It was, therefore, desirable to compare the components of the 
alternating force obtained from the experimental data with those pre- 
dicted by the Hartlen and Currie C3H model. For this purpose the in- 
phase and out-of-phase components of the Hartlen and Currie transverse 
force coefficients have been written equal to the correspond i ng com- 
ponents of the model used herein Csee equation (2)3. This procedure 
yiel ded 

■7— ,A. 2 

327T (^) 

C. (Hartlen and Currie) C . . (29) 

d „ d I 



and 



C (Hartlen and Currie) 
m 



4 ’ 

VT 2 
(— ) 

D 



ml 



(30) 



in which C , and C,, are the drag and inertia coefficients evaluated in 
ml d I 

the present investigation. 

Equations (29) and (30) have been evaluated for A/0 = 0.25 and 0.50 

and the results are plotted as a function of f /f in Figs. 13 and 14. 

c s ^ 
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A similar plot obtained from C3l! is shown in Fig. 15 for A/D = 0.2. 

Even though the A/D values are not exactly identical. Figs. 13, 14 and 

15 could certainly be compared as far as the overall variation of 

and C, are concerned. It is evident that C. curves are fairly similar 
d d 

and that the Hartlen and Currie model could be adjusted sufficiently to 
yield similar C. values for identical A/D values. However, C values 
which are comparable only for < I show opposite variations. It 

has been previously stated that the use of an added mass coefficient 
in terms of the displaced mass of the body is not sufficient and that 
the inertia coefficient of a body oscillating in the transverse direc- 
tion to a uniform flow could be significantly different from its dis- 
placed mass. It suffices to conclude that the Hartlen and Currie model 
does not predict C in accordance with the observations and that the 
overall predictions of the model may be attributed partly to the adjust- 
ment of the other free parameters such as the damping ratio, the density 
ratio, etc. 
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IV. CONCLUSIONS 



The results presented herein have shown that (a) a cylinder can 
undergo self-excited oscillations for certain values of the ratio of 
the natural frequency of the cylinder to the Strouhal frequency of the 
vortices internal damping of the cylinder (?), and the density 

ratio (p^); (b) the drag and inertia coefficients must be determined 
experimentally as a function of A/D and VT/D and must be incorporated 
into the equation of motion; (c) the prediction of the self-excited 
oscillations and their amplitude will be possible with the acquisition 
of additional data for the transverse force; and that (d) the results 
presented herein are quite encouraging and certainly a reflection of 
the combination of the fluid dynamical and vibrational aspects of the 
phenomenon unlike the heuristic models previously proposed. 
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Inertia coefficient versus VT/D for various values of A/D. 





Fig. 4. Drag coefficient versus VT/0 for various values of A/D. 
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Normalized amplitude of the lift force versus VT/D for A/D = 0.50. 
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Normalized amplitude of tt>e lift force versus VT/D for A/D = 0.75, 
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Normalized amplitude of the lift force versus VT/D for A/D = 0.84. 
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Fig. 9. The function I/E versus for various values of the 
damping ratio z,. 
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Fig. 10. Sairiple calculations for the transient and steady-state oscillations of a cylinder 
for the specific parameters indicated. 
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Fig, II, Sample calculations for the transient and sready-state 
oscillations of a cylinder for the specific parameters 
indicated. 
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Fig, 12, Sample calculations for the transient and steady-state 
oscillations of a cylinder for the specific parameters 
i nd icated. 
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Fig. 15. C, and C versus f /f for A/D = 0.25. 
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C. and C versus f /f for A/D 
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Fig. 15. C , and C versus f /f for A/D = 0.2, from 
Hartlen and Currie'" ^C3!]. 
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BFCUTINE RAPINT(X,N,XN,YIM,Y) 



N= NUMBER OF POINTS WHERE THE COORClNdTES CF THE GIVe^ FUNCTION 
ARE KNOWN 

XNtVN = COORDINATES OF THE GIVEN FUNCTION 

X= ABSCISSA OF THE POINT WHOSE ORDINATE WE APE SEEKING 
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1 [F = ID+ 1 

ACX=(X-XN( ID) )/(XN(IDP)-XN( ID)) 
V=VN(IC)+(YN( IDP)-YN(ID) )*AUX 
500 RETURN 
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